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Abstract 

For high-px forward processes at the Large Hadron Colhder (LHC), QCD logarithmic 
corrections in the hard transverse momentum and in the large rapidity interval may both 
be quantitatively significant. The theoretical framework to resum consistently both kinds 
of logarithmic corrections to higher orders in perturbation theory is based on QCD high- 
energy factorization. We present numerical Monte Carlo applications of this method to 
final-state observables associated with production of one forward and one central jet. By 
computing jet correlations in rapidity and azimuth, we analyze the role of corrections to the 
parton-showering chain from large-angle gluon radiation, and discuss this in relationship 
with Monte Carlo results modeling interactions due to multiple parton chains. 



1 Introduction 



Physics in the forward region at hadron colliders is traditionally dominated by soft particle 
production. With the advent of the LHC, forward physics phenomenology turns into a 
largely new field [IH3] involving both soft and hard production processes, because of the 
phase space opening up at high center-of-mass energies. Owing to the unprecedented reach 
in rapidity of the experimental instrumentation, it becomes possible to carry out a program 
of high-pT physics in the forward region. 

Forward jet production enters the LHC physics program in an essential way both for 
QCD studies and for new particle searches, e.g. in vector boson fusion search channels 
for the Higgs boson [HE]. Another area of potential interest in forward physics employs 
near-beam proton taggers [B]: this will enable studies to be made in the central high-px 
production mode with forward protons, which can be used for both standard- candle [7] 
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and discovery physics [S]. In addition to collider physics applications, measurements of 
forward particle production at the LHC will serve as input to the modeling of high-energy 
air showers in cosmic ray experiments [9]. 

The forward production of high-pT particles brings jet physics into a region charac- 
terized by multiple energy scales and asymmetric parton kinematics. In this multi-scale 
region it is compelling to ask [21[in] whether fixed-order next-to-leading calculations reli- 
ably describe the production process or significant contributions arise beyond fixed order 
which call for perturbative QCD resummations. The early observation [11] of potentially 
large logarithmic corrections to jet production at asymptotically high energies has given 
rise to an ample literature of calculations based on the BFKL equation [T214T5] . On the 
other hand, at collider energies both logarithmic corrections in the large rapidity interval 
(of high-energy type) and logarithmic corrections in the hard transverse momentum (of 
coUinear type) are phenomeno logically important. The theoretical framework to resum 
consistently both kinds of logarithmic corrections in QCD perturbation theory is based on 
high-energy factorization at fixed transverse momentum [16]. This factorization program 
is carried through in [T^ for forward jet hadroproduction. 

In this paper we present the application of the results [T7] to the study of jet correlations 
for production of one forward and one central jet at the LHC. The case of forward-backward 
jets will be examined in a forthcoming article. We propose that measurements of hadronic 
final-state observables associated with forward-central correlations can provide, starting 
with the data already taken at the LHC, a sensitive probe of how well QCD multiple 
radiation is taken into account in the Monte Carlo event generators to be used for analyses 
of experimental data in the forward region. The results of such investigations can serve 
to estimate the size of backgrounds from QCD radiation between jets at large rapidity 
separations for Higgs boson searches in vector boson fusion channels. 

Besides the contribution of the higher-order radiative corrections taken into account 
via the results of [T7], the need for realistic Monte Carlo simulations of forward particle 
production raises the question of whether non-negligible effects may come from multiple 
parton interactions. Such multiple interactions are modeled in parton-shower event gen- 
erators used for Monte Carlo simulation of final states at the LHC [T8142T] . and form the 
subject of a number of current efforts [221429] to construct approaches capable of describing 
multiple parton scatterings. In this paper we investigate multi-parton interaction effects 
for forward-central jet correlations within the model [121120] • We observe that, compared 
to the production of multiple jets in the single-scattering picture ^7\, the multi-parton 
mechanism [T8l|20] shifts a significant amount of gluon emissions to larger values of the 
longitudinal momentum fraction x in the initial-state decay chains, because less energy 
is available to each of the sequential parton chains. This results into differences in the 
shapes of the forward-central jet correlations in the azimuth and rapidity plane between 
the single-chain and multiple-chain mechanisms for multi-jet production. This can be in- 
vestigated at the LHC also via measurements of particle and energy flow associated with 
forward production: we leave the study of this to future work. 

The paper is organized as follows. In Sec. [2] we give a concise discussion of the high- 
energy QCD dynamics underlying the hadroproduction of forward jets, based on the results 
of [T7]. We describe the high-energy factorized form of the forward jet cross section that is 
to be coupled to parton showering in order to achieve a full description of the associated 
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hadronic final states. In Sec. Owe discuss aspects of the initial-state parton sfiowers relevant 
to forward hadroproduction. In particular we consider a method to implement parton 
branching at transverse-momentum dependent level not only for gluon-initiated channels 
in the backwards evolution but also for quark-initiated channels. In Sec. |4]we present results 
of numerical Monte Carlo calculations for transverse momentum and pseudorapidity spectra 
and for forward-central correlations in azimuthal angle and pseudorapidity. We compare 
single versus multiple parton interactions, and propose various measurements of forward 
jet observables at the LHC. We give concluding remarks in Sec. [51 

2 Hadroproduction of forward jets 

This section summarizes results from [T7]. In particular we discuss the physical picture 
underlying the factorization formula that will be used for numerical calculations in later 
sections. 

Consider the hadroproduction of a forward jet associated with a hard final state X, 
as depicted in Fig. [T] The kinematics of the process is characterized by the large ratio of 
sub-energies S2/ si ^ 1 and highly asymmetric longitudinal momenta in the partonic initial 
state {xa — !■ 1, a;B — > 0). 




Figure 1: Jet production in the forward rapidity region in hadron-hadron collisions. 

The presence of multiple large-momentum scales in the LHC forward kinematics p|[5D] 
implies that realistic phenomenology of hadronic jet final states requires taking into account 
at higher order both logarithmic corrections in the large rapidity interval (of BFKL type) 
and logarithmic corrections in the hard transverse momentum (of coUinear type). This 
can be achieved via QCD factorization at fixed transverse momentum [16]. A pictorial 
representation of QCD radiative contributions in the rapidity and transverse momentum 
plane is sketched in Fig. |2l Note that kx-factorization is valid to single-logarithmic accuracy. 
In particular, it enables one to obtain logarithmically enhanced terms in rapidity that are 
not associated to any collinear logarithm; conversely, collinear singularities can be taken 
into account to any logarithmic accuracy [21] . This in contrast with calculations in double- 
logarithmic approximations. 

The different expansions in Fig. [2] correspond to different possible ways of reorganizing 
the QCD perturbation series. The results of factorization at fixed kx can be reobtained 
by going to sufficiently sub-leading orders in either the BFKL expansion or the collinear 
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expansion [21]. This applies for instance to the transverse- momentum recoil effects in the 
collinear case [32], and kinematic effects of energy conservation in the BFKL case [33] . 



rapidity 



BFKL 



kj -factorization 



collinear 



Figure 2: QCD radiative contributions to forward jet production in the rapidity and trans- 
verse momentum plane. 



2.1 Factorization of the jet cross section 

Recall that in the case of forward jet leptoproduction [5^ QCD factorization at fixed 
transverse momentum allows one to compute the high-energy asymptotic coefficients for 
the coupling of forward jets to deeply inelastic scattering [16l|35]. Since the early phe- 
nomenological studies [36] forward jet leptoproduction has been investigated at HERA, 
and will play a major role at the proposed future lepton facilities [37] (LHeC, EIC). We 
come back to the possible role of lepton analyses in Sec. O 

In the case of hadroproduction the kx-factorized form of the forward jet cross section is 
given in [17]. This is represented schematically in Fig. [3l Initial-state parton configurations 
contributing to forward jet production are asymmetric, with the parton in the top subgraph 
being probed near the mass shell and large x, while the parton in the bottom subgraph is off- 
shell and small-x. The jet cross section differential in the final-state transverse momentum 
Qt and azimuthal angle (f is given schematically by 

da /" , da ^ 

where ® specifies a convolution in both longitudinal and transverse momenta, a is the hard 
scattering cross section, calculable from a suitable off-shell continuation of perturbative 
matrix elements [T7], ^ajA is the distribution of parton a in hadron A obtained from near- 
collinear shower evolution, and (^g*iB is the gluon unintegrated distribution in hadron B 
obtained from non-collinear, transverse momentum dependent shower evolution. 

Fig. m shows a typical graph contributing to the off-shell matrix element in the qg 
channel. By parameterizing the exchanged momenta k\ and k2 in terms of purely transverse 
four-vectors /cxi and /cx and longitudinal momentum fractions and as 

Pl-P^ = ki= ^ipi + /Cxi + ^lP2 , P2-P6 = k2 = 6P2 + kT+ ^2^1 > (2) 

the forward kinematics implies [IT] (p4 -|- Pq)"^ ^ {ps + p^)'^, ki ~ iiPi, k2 — ^2P2 + ^t, so 
that 

P5-(l-6)Pi , P6 - (1 - 6)P2 - , 6>6- (3) 
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Figure 3: Factorized structure of the cross section. 



In [T7] the full set of the short-distance matrix elements in the forward region, needed for 
the evaluation of the kT-factorized jet cross section ([T]), is computed in exclusive form, for 
all partonic channels. 



Figure 4: A typical graph contributing to the off-shell matrix element in the qg channel. 

As discussed in [17], these matrix elements, though not on shell, are gauge invariant 
and perturbatively calculable. These matrix elements are useful because in the high-energy 
limit they factorize not only in the coUinear emission region but also in the large-angle 
emission region. As long as the factorization is carried out in terms of distributions for 
parton splitting at fixed transverse momentum, they can serve to take into account effects 
of coherence from multi-gluon emission, away from small angles, which become important 
for correlations among jets across long separations in rapidity. We will exploit this in 
performing numerical calculations for central + forward jets in Sec. HI 

2.2 Hard matrix elements and merging 

The precise behavior in transverse momentum resulting from the finite-angle radiation 
taken into account by the method described above is illustrated in Fig. |5] [17]. Here we 
consider the qg channel contribution to forward jet production. (Analogous results for all 
channels can be found in [IT].) We show separately the abelian and non-abelian terms, pro- 
portional respectively to the color factors Cp and CaCf {Cp = {N^ — l)/(2iVc), Ca = Nc). 
With reference to Fig. HI the variables Qt and ip are the final-state transverse momentum 
and azimuthal angle defined by 




P5 




where 
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cos(y9 = Qt ■ ^t/|Qt||^t| , (5) 

while kx is the transverse momentum defined by Eq. ([2]), and effectively measures the 
distribution of the jet system recoiling against the leading di-jets. 




Figure 5: Transverse momentum dependence of the factorizing short- distance matrix ele- 
ments. 

The region kx / Qt — )■ in Fig. [5] corresponds to the leading-order process, with two 
back-to-back jets. The result in this region is simply given by the small-angle limit 



Q'fdd 



^ o^lf'\PT/s) , Qt^Pt= \pt,\ = Iml , (6) 



dQ^dip 

where the function /'•^^ is given by [T7] 

f^'\z) = ^^^==[Clz{l + z)+2C,CA{l~^z + z')] . (7) 

The summation of logarithms for large rapidity y ~ Ins/p^, on the other hand, is 
achieved by convolution of the kx cross section in Fig. [5] with unintegrated parton-splitting 
functions [16]. While in the coUinear approximation case the small- angle result in Eqs. (I6j),(l7l) 
is taken to be valid, in leading order, throughout the range of transverse momentum scales 
up to the factorization scale, we see from the result in Fig. [5] that the role of higher-order, 
multi-gluon emission at large rapidities is to provide a modulation by setting the dynamical 
cut-off at values of kx of order Qt- The essential point is that non- negligible effects may 
arise at high energy from the finite-kx tail. The quantitative evaluation of these effects is 
obtained by integrating the distribution in Fig. [S] over the initial-state parton showers, also 
taken to be transverse-momentum dependent (see Sec. [3D. We will perform such parton 
shower calculations explicitly in the next sections. 

Observe that, as in any parton shower calculation beyond leading order, in order to com- 
bine the hard radiation encoded in the short distance matrix element with the radiation 
from parton showering one needs a specific scheme for merging the two contributions con- 
sistently, by avoiding double counting. The high-energy factorization of Subsec. I2.1l can be 
viewed as providing precisely such a merging scheme. In particular, we recall from 
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that the convolution of off-shell matrix elements over transverse-momentum dependent 
parton-splitting functions is carried out by using systematically the relation for the 2-1-1 
jet cross section 

/ d'kr (^^) ?(M = j d'kr^ - ©(/^ - M ^(Ot)] , (8) 

which provides the necessary small-kx subtraction. 

So in Sees. [3] and H] we will couple Eq. ([1]) to parton showers and perform parton shower 
calculations in the high- energy merging scheme for hadronic final states associated with for- 
ward jets. These calculations will illustrate quantitatively the significance of contributions 
with kj- ^ Qt in the large-?/ region, and will be compared with results of coUinear-shower 
generators, which do not include such finite-kx effects. We note in particular that the 
dependence on the azimuthal angle shown in Fig. [5] is of direct relevance, as we will see 
in Sec. HI for forward-region measurements involving azimuthal plane correlations between 
jets far apart in rapidity. 

It is worth noting that the approach summarized above allows forward jets to be pro- 
duced either from the hard scatter subprocess or from the parton evolution subprocess. We 
will see an explicit numerical illustration of this in Sec. HI This picture can be contrasted 
with the picture from coUinear [10] and BFKL [15] approaches, in which forward jets are 
produced by hard matrix elements or impact factors. This feature of the present approach 
can be traced back to the fact that the factorization ([1]) provides the correct interpolation 
at high energy between the coUinear emission and finite angle regions [TB] . 

Because forward jet production probes the gluon density function for small x (see discus- 
sion around Fig. [3]), it can naturally be used to investigate possible nonlinear effects [39H1T] 
at high parton density. We do not pursue the study of such effects in the present paper; 
but we stress that the formulation [17] at fixed transverse momentum is suitable to de- 
scribe the approach to the high-density region, since, as explained above, it is designed to 
take into account both the effects from BFKL evolution associated with the increase in 
rapidity and also the effects from increasing p^ described by renormalization group, which 
are found to be also quantitatively significant [12H11] for studies of parton saturation. We 
point the reader to e.g. [i5] - [i7] for first Monte Carlo calculations along these lines, and [H] 
for extension to nucleus-nucleus collisions. 

3 Shower evolution and the unintegrated quark den- 
sity 

To obtain a detailed description of the hadronic final states associated with forward jets, we 
need a full parton-shower calculation. This section describes basic features of the showering 
algorithm that we use for such calculations. 

As noted earlier, the factorization formula in Sec. [2] can be used jointly with parton 
showering. Because in the forward kinematics one of the longitudinal momentum fractions 
X in the initial state becomes small, in order to take full account of multi-gluon emission 
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coherence one needs to keep finite- kx terms in the initial-state parton branching |inHS2]0 
We will include these terms according to the CCFM method [5^155] , in particular including 
also showering for quark density channels as explained belowo 

The effect of small-x coherence terms on azimuthal and transverse-momentum jet cor- 
relations and jet multiplicities has been studied in jinUSO], focusing on the case of jet 
lepto-production. It is found that quantitative effects become more significant with de- 
creasing X and decreasing distance in the azimuthal plane between the leading jets. The 
results |l9l[60] suggest that the inclusion of small-x coherence terms in the initial-state 
shower can be relevant in the case of forward jets at the LHCjfl 

It was noted by numerical calculation in [171 130] that for realistic phenomenology of 
forward jets in the LHC kinematics one needs to take into account contributions from both 
quark-density and gluon-density channels. Since CCFM shower evolution has typically only 
included gluon-density terms [55], here we describe how we implement quark channels. For 
the forward jet case that we are interested in, the quark density contributes at fairly large 
values of x. We will thus focus on the valence quark distributional 

We consider the branching evolution equation at the unintegrated, transverse-momentum 
dependent level according to 

zq) 

z)q,q) , (9) 
where q is the evolution scale. The quark splitting function P is given by 

P{z, h) = a, {kl) , (10) 

with (Xs = Cfo^s/'^- Note that, unlike the CCFM kernel given in the appendix B of [51], in 
Eqs. fP]). ffTUl) the non-Sudakov form factor is not included, because we only associate this 
factor to 1/2; terms. The Sudakov form factor is given by 

/ r'}^ dq^ r^~Qo/q i _ \ 
As{qi, Ziqi_i) = exp - / — as(g^(l - zf) dz . (11) 

Here the fractional energy of the exchanged quark i is given by Xi, and the energy transfer 
between the exchanged quarks i — 1 and i is given by Zi = Xi/xi-i. The term xQ^q in 

^This brings in so-called unintegrated, or transverse-momentum dependent, parton distributions. Gen- 
eral issues on these distributions are now actively investigated by many authors, and will influence the use 
of parton branching methods. We comment on this in Sec. [Sj More comments may be found in [49 l l53 j . 

^Alternative methods for taking into account finitc-kx terms in the parton shower are considered in |56l - 
[58] . See [53l|59] for overviews of the subject. 

•^Effects similar to those computed in [49] may also affect azimuthal distributions of h jets [6T1[62] and 
jet multiplicities associated with Higgs boson production j63]. Early measurements of forward jets at the 
LHC can be helpful in this respect to test how well initial state radiation is described and/or for the QCD 
tuning of Monte Carlo event generators. For the counterpart of this in the case of central jets see the first 
LHC measurements [64ll65] . 

"^Work to treat the sea quark distribution at unintegrated level is underway |66| . 



xQy{x,kt,q) = xQvo{x,kt,q) + 



dz 
z 



dq^ n^- 



X As{q,zq)P{z,kt) xQ^ ^-,n.t 



X 



kt + a 
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Eq. on]) is the contribution of the non-resolvable branchings between starting scale go and 
evolution scale q, given by 

xQ^Q{x,kt,q) = xQyQ{x,kt,qo)As{q,qo) , (12) 

where A., is the Sudakov form factor, and the starting distributions at scale go are param- 
eterized using the CTEQ5 u and d valence quark distributions [HZ] as 

xQ^oi^,kt,qo) = xQ^CTEqi^^Qo) exp[-A;^V-^^] > (13) 

with A = 0.92 GeV. 

We next solve Eq. (|9]) numerically. In fig. |6]|^left) the unintegrated u-quark and d-quark 
distributions are shown as a function of x and as a function of kf. In fig. [6l^right) we show 
the following integral of the quark distribution 

xQv{x,kt,q) dkt , (14) 

and compare this for consistency with the distribution obtained from CTEQ [67] at the 
same scale. 




Figure 6: Left: Unintegrated quark distribution (u,d-quarks) as a function of x at fixed kt 
(top) and as a function of kt at fixed x (bottom) at a scale g = 4 GeV. Right: Integral of 
the unintegrated quark distribution (u,d-quarks) as a function of x for different scales q. 
Also shown is the u, d-quark distribution obtained from CTEQ fU^ . 



9 



4 Central + forward jet production at the LHC 

In a typical LHC experiment jets can be measured for high transverse energy > 30 GeV 
in a large range of pseudorapidity 77. In the following we consider differential cross sections 
for dijets (Fig. [7]) reconstructed with the Siscone algorithm pH] with R = 0.4, where one 
jet is in the central region defined by |?7c| < 2 and the other jet is in the forward region 
defined by 3 < |?7/| < 5. 



away from 
the jets 




between the 
jets 



forward 
jet 



Figure 7: ( top ) Jets in the forward and central detectors, and azimuthal plane segmentation; 
(bottom) particle and energy flow in the inter-jet and outside regions. 



4.1 Dijet cross section and /S.R distribution 

The total cross section for a central and a forward jet obtained with the Cascade [SI] 
Monte Carlo event generator (version 2.2.03 including the matrix element calculated in |17j . 
the unintegrated gluon distribution set A and the unintegrated valence quark distribution 
described in section [3]) is given in Tab. [H We compare the prediction from Cascade with 
the prediction from the Pythia [20] Monte Carlo event generator running in two modes: 
with and without multi-parton interactions (Fig. [8]). We use tune PI [18], which allows for 
more radiation from parton shower. Both Monte Carlo generators simulate higher order 
QCD corrections with parton showers: CASCADE uses parton showers according to the 
CCFM evolution equation whereas Pythia uses DGLAP based parton showers. The total 
cross section predicted by Cascade lies in between the prediction of Pythia with and 
without multiparton interactions (Tab. [1]). 

In Fig. [H] we investigate the physical mechanism producing the central and forward 
jets. We plot Ai? = ^(A0)2 + (A?7)2, where A0 = (jy^et - 4>part (A?7 = r]jet - Vpart) is 
the azimuthal (pseudorapidity) difference between the jet and the corresponding parton 
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Figure 8: Multi-jet production by (left) multiple parton chains; (right) single parton chain. 

Table 1: Integrated dijet cross section for > 10(30) GeV in the range |?7c| < 2 and 
3<|^/|<5. 





a{ET > 10 GeV) 


a{ET > 30 GeV) 


Cascade 


469 /ib 


3.1 /ib 


Pythia (MPI) (PI) 


798 /ib 


3.5 /ib 


Pythia (no MPI) 


346 /ib 


3.3 /ib 



from the matrix element. The AR distribution allows one to see whether the jets are 
dominated by hard partons from the matrix element, or whether they originate from the 
parton shower. For the low Et jets the distribution in AR has a significant contribution 
from jets not corresponding to a parton from the matrix element {AR > 1). 

The bump structure of the distribution at AR > 1 is consistent with random distribu- 
tions in A(f) and Arj within the phase space region investigated here, showing that indeed 
AR > 1 corresponds to the region where there is no correlation between the parton from 
the matrix element and the jet. 

It is interesting to observe that CASCADE predicts a similar distribution as obtained 
from Pythia with multiparton interactions, whereas Pythia without multiparton inter- 
actions has a significantly smaller contribution at large AR. The situation changes for the 
high Et jets: CASCADE predicts significantly more jets (especially in the central region) 
coming from the parton shower as compared to Pythia. This behavior is understandable, 
as the small-x initial-state parton shower allows for higher transverse momentum radia- 
tion compared to a collinear parton shower. It also shows, that the most forward jet is 
mainly coming from the matrix element parton, whereas the central jet has a significant 
contribution from the parton shower. 

4.2 Transverse momentum spectra 

Figs. [To] and [TT] show the differential cross section da/dEx for jets reconstructed with the 
Siscone algorithm in the central region defined by |?7c| < 2 and in the forward region defined 
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Figure 9: AR distribution of the central (\r]c\ < 2, left) and forward jets ( 3 < \7]f \ < 5, 
right) for Et > 10 GeV (upper row) and Et > 30 GeV (lower row). The prediction from 
the k± shower ('CASCADE^ is shown with the solid blue line; the prediction from the collinear 
shower (T*YTHIA ) including multiple interactions and without multiple interactions is shown 
with the red and purple lines. 
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by 3 < |?7/| < 5. The left(right) plots show the results when both jets have > 10(30) GeV. 
We see that the k^-dependent parton shower implemented in Cascade produces a signif- 
icantly harder spectrum especially for the jet in the central region when both jets are 
required to have Et > 10 GeV. The predictions from Pythia with and without multipar- 
ton interactions are similar at large Erp, whereas the multiparton interactions contribute 
significantly in the low Ej- region {Et < 30 GeV) resulting in a larger cross section. 




50 100 150 200 

Et (GeV) 




50 100 150 200 

Et (GeV) 



Figure 10: Transverse momentum spectra of central jets with \ric\ < 2 at y/s = 7 TeV for 
Et > 10 GeV (left) and for Et >30 GeV (right) for events which have a forward jet with 
Et > 10(30) GeV in 3 < \r]f\ < 5. The prediction from the kj_ shower ('CASCADEj is 
shown with the solid blue line; the prediction from the collinear shower (T'YTHIAj including 
multiple interactions and without multiple interactions is shown with the red and purple 
lines. 



This behavior can be understood since CASCADE uses matrix elements which are calcu- 
lated within high-energy factorization, allowing harder transverse momentum dependence 
as compared to collinear factorization. Moreover Cascade uses the CCFM parton shower 
with angular ordering which at small x allows for a random walk in transverse momentum, 
and thus allows for more and harder parton radiation compared to DGLAP based parton 
shower as implemented in Pythia. 

It is also interesting to observe that the contribution from multiparton interactions as 
implemented in Pythia is important only for the region of low Et jets. We have checked 
also the prediction obtained from Pythia using a different parameter set for the modeling 
of multiple-parton chains (tune D6T [I9]). We find that the difference between the different 
tunes for the multiparton interaction parameters is smaller than the difference coming from 
the noncoUinear corrections to single-chain parton shower. 

To see better the significance of the results in Figs.[TO]and[TTl we have also considered the 
next-to-leading-order Monte Carlo generator Powheg [69]. We find that the Et spectra 
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Figure 11: Transverse momentum spectra of forward jets with 3 < |?7/| < 5 at y/s = 7 TeV 
for Et > 10 GeV (left) and for Et > 30 GeV (right) for events which have a forward 
jet with p_L > 10(30) GeV in \r]c\ < 2. The prediction from the k± shower ('Cascade^ is 
shown with the solid blue line; the prediction from the collinear shower fPYTHlA^ including 
multiple interactions and without multiple interactions is shown with the red and purple 
lines. 
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from PowHEG, in the range of rapidity and transverse energy considered here, are very 
close to those of Pythia, indicating that the enhancement from the shower in Figs. [10] 
and [11] is not simply due to the next-to-leading term, but comes from corrections beyond 
next-to-leading order. This feature is noteworthy; its study should be extended to a wider 
range of forward-jet observables. 

4.3 Rapidity dependence 
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Figure 12: Pseudorapidity spectra of produced jets for ^/s = 7TeV with requirement that 
Pt > 10 GeV (left) andpT>30 GeV (right). Cascade is shown with the blue solid line, 
Pythia with (without) multiparton interactions is shown as the dashed red (dotted purple) 
line. 

In Fig. [T2]we show the differential cross section da/dt] for dijet events with Et > 10(30) 
GeV in two regions of < |?7| < 2 and 3 < \ri\ < 5. Again we compare the prediction from 
Cascade with the one from Pythia without and with multiparton interactions. 

We observe that the cross section in the central region for Cascade is rising towards 
larger rj whereas for Pythia the cross section is flat. The cross section in the forward region 
is steeply falling towards large 77. The slope of this distribution is different from Cascade 
and Pythia. A closer investigation of these different behaviors is underway. A significant 
contribution to the difference comes from the treatment of the quark distribution (Sec. 3), 
and suggests the need to include both valence and sea quark distributions at unintegrated 
level. 

From Fig. [12] (left) we see that at small Et the influence from multiparton interactions 
in Pythia is significant. Cascade predicts a cross section of similar size as Pythia with 
multiparton interactions in the region 3 < 77 < 3.5, but the distribution falls more rapidly 
towards larger rj. In the large Et region (Fig. [12] (right)) the differences in cross sections 
as a function of rj become smaller. 
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4.4 Azimuthal dependence 

The azimuthal correlation of a central and forward jet is a measure of the parton radia- 
tion between the jets and is therefore a probe of how well the Monte Carlo parton shower 
is simulating the higher-order parton emissions. The azimuthal decorrelation of forward 
and backward jets has been proposed as one of the measurements to test BFKL dynam- 
ics [in] (and references therein). In lowest order, BFKL predicts a much larger decorrelation 
compared to calculations in coUinear factorisation. With a larger separation of the jets in 
A?7, the phase space for parton radiation is increased. However, significant multiparton 
interactions could perhaps mimic a signal expected from small x dynamics. 

In Fig. [13] we show the differential cross section d^cr/dA^A?]. The decorrelation as a 
function of Arj increases in Cascade as well as in Pythia. In the low Ej^ region (Fig. 1131 
(left)) the increase in decorrelation with increasing At] is very significant. The cross sec- 
tion for jet separation up to At] < 4 is very similar between Cascade and Pythia with 
multiparton interactions, whereas a clear difference is seen to Pythia without multipar- 
ton interactions. However, at large At] > 4 the decorrelation predicted by Cascade is 
significantly larger than the prediction including multiparton interactions. 

In the higher Et region CASCADE predicts everywhere a larger decorrelation. In this 
region, the influence of multiparton interactions in Pythia is small and the difference to 
Cascade comes entirely from the different parton shower. 

In Fig [m we show the average (cos(A0 — tt)) as a function of the rapidity separation 
Arj of the central and forward jets. This quantity is considered in [TU |[T3] - [T3] as a sensitive 
probe for BFKL dynamics. We observe that the distribution of (cos(A0 — vr)) at pj- > 
10 GeV shows only little difference between Cascade and Pythia, whereas the differential 
distribution of dPa / dAipArj (Fig. [T3|) is more discriminative. In the px > 30 GeV case 
Cascade predicts a larger decorrelation than Pythia, consistent with what is observed 
in Fig. M 

5 Summary and outlook 

The production of hadronic jets in the forward region of pp collisions will form a largely new 
area of experimental and theoretical activity at the Large Hadron Collider. Forward jet 
production will enter the LHC physics program both for new particle discovery processes 
(e.g., vector boson fusion channels for Higgs boson searches) and for new aspects of standard 
model physics (e.g., QCD at small x and its interplay with cosmic ray physics). 

In this paper we have focused on the study of jet correlations for production of a forward 
and a central jet at the LHC. The capabilities of LHC forward and central detectors allow 
measurements of such correlations to be made for the first time across large rapidity inter- 
vals {Ay > 4 6). In this kinematic region the evaluation of QCD theoretical predictions 
is made complex by the presence of multiple mass scales, which raises questions on the 
potential need for, on one hand, perturbative QCD resummations and, on the other hand, 
possible corrections beyond single parton scattering. 

To address the structure of the final states associated with forward jet production, 
in this paper we have used a merging scheme based on high energy factorization [T6l[T7] 
to combine hard matrix elements and parton showering. This is designed to incorporate 
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Figure 13: Cross section as a function of the azimuthal difference between the central 
and the forward jet for different separations in pseudorapidity Ar], at y/s = 7TeV for jets 
with pt>10 GeV (upper) andpT>SO GeV (lower). Cascade is shown with the blue solid 
line, Pythia with (without) multiparton interactions is shown as the dashed red (dotted 
purple) line. yj 
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Figure 14: Average (cos(A(/) — vr)) between the central and the forward jet as a function 
of separation in pseudorapidity Arj, at y/s = 7TeV for jets with px > 10 GeV (left) and 
> 30 GeV (right). Cascade is shown with the blue solid line, Pythia with (without) 
multiparton interactions is shown as the dashed red (dotted purple) line. 



radiative contributions in such a manner that both logarithmic corrections in the large 
rapidity interval and logarithmic corrections in the hard transverse momentum are taken 
into account to higher orders in Ug- Both kinds of contributions are likely to be needed 
for reliable phenomenology of forward hard production. Their summation is achieved by 
including finite-kx terms in the parton branching and matrix elements. 

In this approach forward jets may be produced from either the hard scatter subprocess or 
the parton showering subprocess. The resulting physical picture of forward hard production 
is thus different from that of purely BFKL or collinear calculations, in which, also at the 
next-to- leading order (see e.g. [IOl[T5]), forward jets are produced by hard impact factors 
or matrix elements. 

We have implemented the high-energy merging scheme in the hadron-level shower Monte 
Carlo event generator Cascade [M], and used this to obtain numerical predictions for 
several forward-jet observables. In doing this, we have also implemented in the Monte 
Carlo generator [51] an algorithm to include the parton branching evolution for the valence 
quark distribution at unintegrated level, which extends the previous CCFM algorithm [55] 
for the unintegrated gluon distribution. 

We have examined the effects of the higher-order radiative contributions taken into 
account by this approach by computing the jet transverse-momentum and rapidity spectra 
and the jet correlations in rapidity and azimuth. We find that the effects are significant 
especially in the slope of the Et spectrum and in the jet angular correlations. In particular, 
we find that while the average cosine of the azimuthal separation between the leading jets 
is not affected very much as a function of rapidity by finite-kx terms, the detailed shape of 
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the A0 distribution is. 

We have also investigated the model [18], [20] for multiple parton interactions, corre- 
sponding to corrections beyond single parton scattering. Our analysis shows that certain 
features of forward jet production such as rapidity and Et spectra, found by including 
high-energy, noncoUinear corrections to single-chain parton showers, can be mimicked by 
effects of multiple parton chains. However, distinctive shapes are found both in the AR 
distribution and in the azimuthal correlations. We suggest that measurements of the parti- 
cle and energy flow in the region both between the jets and away from the jets should have 
stronger discriminating power between the single-chain and multiple-chain mechanisms for 
multi-jet production. The detailed analysis of this point will be the subject of a sepa- 
rate paper. Note that this also points to the phenomenological relevance of energy flow 
observables such as those investigated in [TOHT^] . 

As observed in Sees. [2] and [3l many of the theoretical tools that underlie forward jet 
physics, from parton branching beyond leading order to perturbative QCD resummations 
to, possibly, the approach to the saturation region, depend on the notion of transverse 
momentum dependent, or unintegrated, parton distribution functions (u-pdfs). In the 
calculations of this paper we take the high-energy definition of u-pdfs [TB], namely, we rely 
on the fact that for small x u-pdfs can be defined gauge-invariantly (and can be related 
to the ordinary pdfs renormalized in the minimal subtraction scheme MS [73]) by going to 
the high-energy pole in physical amplitudes [16]. More general characterizations of u-pdfs, 
valid over the whole phase space, are very desirable. A complete framework is yet to be 
fully developed though. Recent results in this area, see e.g. [71H5D] . are likely to eventually 
have a bearing on forward jet physics. 

Note that forward hard production processes will be relevant not only for LHC physics 
but also for physics at the planned future lepton facilities [37] (LHeC, EIC). Thus a unified 
understanding in hadron-hadron and lepton-hadron collisions is desirable. As recalled in 
Sec. 121 QCD high-energy factorization [16] has been used to determine the asymptotic 
coefficients [35] that couple forward jets to deeply inelastic scattering. Since the early 
phenomenological studies [36], forward jet leptoproduction has been investigated at HERA. 
Measurements of forward jet cross sections at HERA [81] have illustrated that neither 
fixed-order next-to-leading calculations nor standard shower Monte Carlo generators [T2| 
IHI1IS2], e.g. Pythia or Herwig, are able to describe forward jet ep data. This provides 
additional motivation for developing methods capable of treating the multi-scale kinematics 
and describing jet production beyond the central rapidity region. It is of interest to analyze 
HERA data [83] looking at forward + central jets, similarly to what is done in this paper 
(Sec. H]) for pp collisions; but due to the phase space available for multiple jet radiation, 
such studies are likely to prove much more relevant at a future high-energy lepton collider. 

The analysis performed in this paper shows that the final states associated with forward 
jet production at the LHC receives significant contributions from radiative corrections that 
take into account both large logarithms of rapidity and large logarithms of transverse 
momentum. Distinctive effects are found, for instance, for the distribution in the azimuthal 
separation A0 between forward and central jets. We obtain distinctive predictions both 
with respect to parton showers modeling multiple parton interactions and with respect 
to parton showers including next-to-leading fixed-order corrections. This analysis can be 
extended to correlations of forward and backward jets. It can thus serve to estimate the 
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size of backgrounds from QCD radiation in new particle searches from vector boson fusion 
channels. 
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